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We consider the proximity effect in a normal dot coupled to a bulk superconducting reservoir by 
the tunnel contact with large normal conductance. Coulomb interaction in the dot suppresses the 
proximity minigap induced in the normal part of the system. We find exact expressions for the 
thermodynamic and tunneling minigaps as functions of the junction's capacitance. The tunneling 
minigap interpolates between its proximity-induced value in the regime of weak Coulomb interaction 
to the Coulomb gap in the regime of strong interaction. In the intermediate case a non-universal 
two-step structure of the tunneling density of states is predicted. The charge quantization in the 
dot is also studied. 
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When a normal metal forms a contact with a supercon- 
ductor it acquires some superconductor features. One of 
them is the suppression of the electron density of states 
(DoS) at the Fermi energy. This effect is governed by the 
Andreev processes at the normal-metal-superconductor 
(NS) boundary. Cooper pairs tunnel into the normal 
metal preserving their phase and thus inducing the su- 
perconductive correlations. Coulomb repulsion leads to 
phase fluctuations and reduces this proximity effect 0. 
For a disordered normal metal, a minigap of the order of 
the inverse escape time appears in the excitation spec- 
trum . Below we study the suppression of the minigap 
by Coulomb interaction. 

To measure the density of states the tunneling spec- 
troscopy technique is widely employed. The conductance 
as a function of the voltage between the external probe 
and the system under investigation is proportional to the 
tunneling density of states (TDoS). Coulomb interaction 
affects the TDoS suppressing tunneling conductance at 
small voltages. This suppression is known as the zero 
bias anomaly 0, 0| . The interplay between the zero bias 
anomaly and the proximity effect was first studied in 
Ref. 5] for a 2D thin normal film coupled by the tun- 
nel junction to a superconductor. The renornialization 
group procedure yields the power-law suppression of the 
minigap in the TDoS by the Coulomb interaction. 

Coulomb repulsion is ultimately strong in restricted 
geometry. In such systems direct observation of charge 
quantization is made possible by the Coulomb blockade 
effect 0. When a zero-dimensional (OD) normal grain is 
coupled to a normal reservoir by the tunnel junction with 
small dimensionless [in units of e^/?i] conductance G, the 
equilibrium charge of the grain is a step function of the 
gate voltage at zero temperature. Contrary, if G » 1 
the charge of the grain is no longer conserved and only 
exponentially weak modulation of charge- voltage depen- 
dence remains 0- The situation changes dramatically if 



the reservoir is superconducting . Now a single electron 
cannot escape into the superconductor and charge quan- 
tization is observed even at large normal conductance. 

In this Letter we consider the OD NS system described 
above. The parameters of the system are assumed to 
satisfy the following conditions: 



£;Th> A»(£;g,i;c)»^. 



(1) 



Here .ETh is the Thouless energy, A is the superconductor 
gap, 5 is the electron mean level spacing in the grain per 
one spin component. Eg — GS/A is the proximity mini- 
gap 2] in the absence of Coulomb repulsion, and Ec = 
e^/(2G) is the charging energy, with G being the junc- 
tion's capacitance. The dimensionless conductance G is 
assumed to be large. We neglect effects of quasiparticle 
transport assuming temperature sufficiently low and thus 
quasiparticle conductance Gqp = Gexp(— A/T) ^ 1. 

Coulomb interaction modifies the thermodynamic and 
tunneling DoS in different ways. The minigap Eg in the 
thermodynamic DoS is gradually suppressed with the in- 
crease of Ec due to enhanced phase fluctuations, whereas 
the dependence of the minigap E'*"" in the TDoS on Ec 
is more complicated. Being suppressed at weak inter- 
action by phase fluctuations, it eventually reaches the 
Coulomb gap Ec in the strong interaction limit. Qual- 
itatively, these two regimes are distinguished by the re- 
lation between the charging energy Ec and the energy 
Ej = {Eg/S)\og{A/Eg), which is the Josephson energy 
of the fictitious system where the normal grain is replaced 
by the weak superconductor with the gap Eg. Below we 
find the exact dependence of both the DoS and TDoS on 
the strength of the interaction for arbitrary Ec/Ej. 

We also study the charge quantization in the grain ac- 
counting for both Coulomb and proximity effects. The 
result of Ref. |1] is valid only for Ec Ej. In the 
opposite case proximity coupling smears the Coulomb 
staircase preserving however the sharp discontinuities at 



2 



half-integer charge. 

To investigate the electron properties of the NS sys- 
tem we use the non-linear cr-model in the Matsubara 
representation. Disorder average is made with the help 
of replica trick. Thus the cr-model is formulated in terms 
of matrix field Q operating in the Nambu-Gor'kov (Pauli 
matrices Ti), energy, and replica spaces. The OD Coulomb 
interaction is decoupled by the electric potential (f). The 
cr-model action reads 01 



^[Q, 0] = - J Tr [(era -f 0)Q]-^ Tr[Q50]- 



dr- 



(2) 

Hereafter we omit replica indices for brevity. We will use 
fi for the Q-matrix of the superconductor, Qs, restrict- 
ing our model to the subgap re gion £ ^ A. The contribu- 
tion from higher energies leads to the renormalization 
of the capacitance C > C + e'^G/ (2A) and corresponding 
renormalization of the charging energy Ec- 

In the action (O the matrix Q is linearly coupled to 
the potential (j). This means that Q contains not only 
soft electron modes of the system with energies close to 
the Fermi energy but also the high energy fluctuations 
corresponding to the shift of the whole energy band by 
electric potential. To get rid of this contribution we make 
a gauge transformation proposed in Ref. [lH: Qrr' — 

^^T^K{r)Q^^^^~^HK{r') ^j^j^ ^(^^ ^ f (l,{T)dT. In tCrmS 

of these new variables the action has the form 



S[Q,K] 



■Tr(ef3Q) 



dr 



2ttE„ 



AE, 



c 



Qi^} cos2K + Q^^'J sm2K 



(2), 



(3) 



where the symbols Q*-*-* = tr(fi(5)/2 denote the Nambu- 
Gor'kov components of the matrix Q. 

The general dynamics governed by the action Q is 
complicated as the variables Q and K are strongly cou- 
pled with each other. Fortunately, in the region of the 
parameters specified by Eq. ^ the characteristic fre- 
quencies of the variable Q (of the order of the renor- 
malized proximity gap Eg) appear to be much smaller 
than those of the variable K. Therefore, it is possible 
to employ the adiabatic approximation, integrating first 
over the "fast" variable K. The resulting effective action 
S[Q] for the "slow" variable Q should then be subject to 
the saddle-point approximation (SPA). Fluctuations of 
Q-matrix describe the effects of level statistics, therefore 
the SPA works fine if the relevant energy scale Eg exceeds 
6. Such a condition is definitely satisfied in the regime 
of weak Coulomb repulsion: Eg/S = G/4 ^ 1. Below 
we will find the applicability range both of the adiabatic 
and saddle-point approximations at arbitrary strength of 
the Coulomb interaction. 

The action averaged over K is obviously uniform in 



time that allows to assume the stationary saddle point 



Qe 



2TrS{e — e'){T3 cosa(e) -I- fi sina(e)), (4) 



parametrized by the single function a{e). Possibility to 
choose Q in the form without the T2 term is related 
to the zero mode with respect to K i-^ K + const. 

Substituting the ansatz into the action ^ we find 
the imaginary time evolution for K. It is determined by 
the simple Hamiltonian 



H ^Ec [-d^/dK'^ - 2gcos2if] 



(5) 



The parameter q ~ {Eg/2EcS) J_^desma{e) measur- 
ing the relative strength of the proximity coupling with 
respect to Coulomb interaction is to be determined self- 
consistently later on. The 27r-periodic boundary condi- 
tions for K are implied. 

At zero temperature the dynamics of K(t) is frozen 
at the ground state of 10. The ground state energy can 
be expressed in terms of zeroth Mathieu characteristic 
value aQ{q): Eq = Ecao{q). After averaging over K one 
obtains the effective action for Q: 



S[Q] 



J dr — i J de e cosa{e) + Eo{q) 



(6) 



Minimizing this expression we find the BCS-like solution 
a{e) = arctan(£'g/e). The renormalized value of the 
minigap Eg satisfies the following system of equations: 



1 dEp 
2Ec dq ' 



Q = 



EgEg 2A 



(7) 



The dependence of the ground state energy Eq on q is 
especially simple in two limiting cases. If the Coulomb 
interaction is weak the parameter q is large and we can 
estimate the potential energy in Eq. (O by the oscilla- 
tor potential. This gives a small correction to the bare 
minigap value Eg: 



En — E„ 



EcS 



\og{2 A/ Eg)' 



g > 1. 



(8) 



In this regime the variable K has the frequency Ec^/q 
which is yjEc/5 times larger than Eg. In the op- 
posite limit of strong Coulomb interaction the poten- 
tial in Eq. jSJ can be treated perturbatively yielding 
00(9) — "9^/2- The minigap appears to be exponen- 
tially small in Ec'- 



Eg = 2A exp 



2Ec5 
El 



g < 1. 



(9) 



At the same time, the frequency of K is Eq ^ Eg, en- 
suring the validity of the adiabatic approximation. 
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FIG. 1: The minigap in the thermodynamic DoS as a func- 
tion of the Coulomb energy Ec- Dashed hnes illustrate the 
asymptotic behaviors Q and Q. The crossover from weak to 

strong fluctuations regime occurs at Bc^/Sg ~ log(A/S9)/2. plotted assuming log(A/£;g) = 5 and G = 40. Inset: two- 



FIG. 2: The minigap in the tunneling DoS as a function of 
the Coulomb energy. For small Ec it follows Eg given by 
Eq. Q. In the opposite limit S*"" = Ec- This figure 



This figure shows the case log(A/£'g) = 5. 

The whole dependence of Eg on Ec is shown in Fig. ^ 
The thermodynamic DoS itself is just the conventional 
BCS density of states with the gap Eg. 

Now we turn to the tunneling density of states. In the 
cr-model formalism the TDoS is the analytic continuation 
of the function 



1 



p{e) = - / dre'^Ur(f3e'"^'"ge- 

J \ I Tfl 



(10) 



from positive Matsubara energies e to real energies E. 
The angular brackets imply averaging with the action ^ . 
In the adiabatic approximation this means averaging 
over K with the Hamiltonian ||SJ) and substitution Q 
Q. The result for the TDoS is the convolution of the 
thermodynamic DoS with the phase correlator C(t) = 
{e'iK{r)~Km^ given by 



Ciiu) 



E 



2{En ~ Eq) 
{En - Eof + 



(11) 



in the frequency representation at zero temperature. 
Here Ei are the energy levels of the Hamiltonian jSJ. 
The matrix elements in Eq. (|ll|l are non-zero only for 
n ~ Ak + 1 and n — 4fc + 2. Evaluating the convolution 
of the thermodynamic DoS with C{lu) and performing 
analytic continuation, we obtain for the TDoS: 



piE) = 



|2 |x„| 0{\Xn\ - Eg 



El 



(12) 



where Xn = E — (£'„ ^ Eq). The tunneling minigap as- 
sociated with the n — I term is given by 



T^tun 



Er, 



El — Eq. 



(13) 



Fig. 121 shows its dependence on Ec which interpolates 
from the proximity gap (jHJ in the weak interaction regime 



step tunneling density of states in the intermediate regime. 
The charging energy is such that Ec5/Eg — 2.8 (as shown by 
the dot on the main graph) which corresponds to g = 0.95. 



to the Coulomb gap Ec at large Ec. Note the reentrant 
behavior of Eg^'^ as a function of Ec at small Ec. 

The form of the TDoS depends on the strength of the 
Coulomb interaction. Due to the fast decrease of the 
matrix element |(0|e*^|n)| with n only the first two 
terms may significantly contribute to the sum in Eq. I|12|l . 
In the weak Coulomb interaction regime (g ^ 1), the 
n — I contribution dominates. Since the level split- 
ting El — Eq is exponentially small the TDoS coincides 
with the thermodynamic DoS in this limit. When the 
Coulomb interaction is strong (q <^ 1), the matrix el- 



1) 



(0|e^^|2) 



1/2 and p(E) 



ements KUIe" 

{2/S)d{\E\ — Ec). In the intermediate regime (q ^ 1) 
the TDoS acquires a non-universal two-step structure, 
being zero for E < i?^"", of the order of [(Ole'-^ll)]^ 
for £;^"" < E < Eg + E2 - Eq, and close to 1 for 
E > Eg + E2~EQ, see the inset of Fig. [3 The BCS-type 
singularities are still present at the step edges but they 
are relatively weak. 

At non-zero temperature all above results are still valid 
provided T <^ Eg. If the temperature is higher all inte- 
grations over energies should be replaced by summations 
over Matsubara frequencies. Moreover, the phase K is 
not frozen at the ground state in this case, and the tem- 
perature dependent free energy should be used instead of 
the ground state energy in Eq. ^ and below. These cal- 
culations show that the true tunneling minigap vanishes 
when the temperature is of the order of Eg{T — 0). 

Now we briefly comment on the validity of our ap- 
proach. The above results were obtained using the adi- 
abatic approximation for the integral over K and the 
saddle point approximation for Q. To justify this ap- 
proach one can calculate the first fluctuation correction 



4 



o 



0.5 

Eg /Eg 



20 




Ec{n — NY . The second order in q correction to the 



FIG. 3: The dependence of thermodynamic minigap, tunnel- 
ing minigap, and average charge on the gate voltage V (ver- 
tical axis) at three different values of Coulomb energy. The 
parameters are G — 40, log(A/i5g) — 5. 



to the TDoS. To this end we parametrize Q in the stan- 
dard way in terms of the matrix W, and average over 
K with the Hamiltonian (0) taking into account W per- 
turbatively. As a result we end up with the quadratic 
action for the fluctuating W. Evaluating the first order 
perturbative correction to the density of states we find 
that it is small provided that Eq ^ S (adiabaticity con- 
dition) and Eg S (validity of the SPA). All the above 
results were obtained in the lowest-order approximation 
over tunneling probability T across NS interface. The 
effect of higher-order terms can be neglected if Andreev 
conductance through this interface is small, Ga <C 1. 

Finally, we consider the case of a non-zero gate voltage 
V applied to the normal grain in the standard Coulomb 
blockade setup 0. If the grain is coupled to the gate 
by the capacitance Cg the charging energy becomes 
Ec{Q/e — N)^, where Q is the charge operator of the 
grain, eiV = CgV is the equilibrium charge, and we re- 
defined Ec = eV2(C + Cg). The Hamihonian for K is 
changed to 



H ^ Er 



-id/dK- N) - 2gcos2X 



(14) 



The energy spectrum of this Hamiltonian has the stan- 
dard band structure. The boundary conditions pick out 
two energy levels from each band with the quasimomen- 
tum determined by the gate potential. When the equilib- 
rium charge N approaches half-integer values these two 
energy levels cross. Since all physical quantities depend 
periodically on N, in what follows we assume |7V| < 1/2. 

In the weak Coulomb interaction regime, Ec ^ Ej, 
the lowest energy band is exponentially narrow. The 
lowest level of the Hamiltonian depends weakly on N, 
and all the results for the DoS and TDoS obtained 
above are left unchanged. For strong Coulomb inter- 
action we calculate the ground state energy perturba- 



ground state energy Eq 



(2) 



-Ecq'^/2{1 - N'^). From 



the self-consistency equations Q we find the thermody- 
namic minigap Eg = 2Aexp[-{2EcS/ E'^){1 ~ N'^)]. The 
minigap in the TDoS is determined by Eq. H13I) yield- 
ing El""^ ^ Eg + Ecil-2\N\). This quantity depends 
strongly on TV decreasing to the exponentially small value 
Eg at N = 1/2. These results are illustrated in Fig.|3| 

The average charge of the grain is Q = eN — dEo/dV. 
In the limit of strong Coulomb interaction the re- 
sult of Q is reproduced. In the opposite limit the 
Coulomb staircase is smeared up to an exponentially 
weak modulation. Employing the tight binding approx- 
imation for the Hamiltonian (|14|l we get Q/e = N — 
Cg/(C-H Cg)8V2^g3/4g-4yg sin(7riV). The smah steps of 

the height Cg/{C + Cg)16V2^q^^'^e-'^-^ at half-integer 
charge are still present as predicted in Ref. This 
feature is the consequence of the double-electron charge 
transport through the junction. The Coulomb staircase 
is shown in the right part of Fig. |3 

In conclusion, we have shown that interplay between 
proximity and charging effects in a SIN structure with 
large normal conductance G is governed by the ratio of 
the charging energy Ec to the fictitious Josepshson en- 
ergy Ej = G^S log(A/G(5). DoS and TDoS are calculated 
as functions of Ec/Ej and the gate voltage V. An un- 
sual two-step shape of the TDoS is found at Ec ^ Ej. 
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tively in small q. Neglecting the q term we have En 



(0) 
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